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ABSTRACT 



Context. Low-energy cosmic rays are the dominant source of ionization for molecular cloud cores. The ionization fraction, in turn, 
controls the coupling of the magnetic field to the gas and hence the dynamical evolution of the cores. 

Aims. The purpose of this work is to compute the attenuation of the cosmic-ray flux rate in a cloud core taking into account magnetic 
focusing, magnetic mirroring, and all relevant energy loss processes. 

Methods. We adopt a standard cloud model characterized by a mass-to-flux ratio supercritical by a factor of ~ 2 to describe the 
density and magnetic field distribution of a low-mass starless core, and we follow the propagation of cosmic rays through the core 
along flux tubes enclosing difl'erent amount of mass. We then extend our analysis to cores with difl'erent mass-to-flux ratios. 
Results. We find that mirroring always dominates over focusing, implying a reduction of the cosmic-ray ionization rate by a factor 
of ~ 2-3 over most of a solar-mass core with respect to the value in the intercloud medium outside the core. For flux tubes enclosing 
larger masses the reduction factor is smaller, since the field becomes increasingly uniform at larger radii and lower densities. We also 
find that the cosmic-ray ionization rate is further reduced in clouds with stronger magnetic field, e.g. by a factor ~ 4 for a marginally 
critical cloud. 

Conclusions. The magnetic field threading molecular cloud cores aff'ects the penetration of low-energy cosmic rays and reduces the 
ionization rate by a factor 3-4 depending on the position inside the core and the magnetization of the core. 

Key words. ISM:cosmic rays -ISM: clouds, magnetic fields 



1. Introduction 

Cosmic rays (CRs) are a primary source of ionization in inter- 
stellar clouds (Hayakawa et al. 1961, Spitzer & Tomasko 1968), 
contribute to the heating of the gas (Glassgold & Langer 1973) 
and affect the chemistry of dense clouds by removing molecules 
from the most volatiles ices like CO (Hasegawa & Herbst 1993). 
For a review of CR induced chemistry in the interstellar medium 
(ISM), see Dalgarno (2006). In addition, the observed diffuse 
gamma-ray emission from the galactic plane is believed to be 
the result of the decay of neutral pions produced during inelastic 
collisions of high-energy (> 1 GeV) CRs with the ISM (see e.g. 
Gabici et al. 2007). 

CR protons and heavy nuclei in the range 1 MeV/nucleon- 
1 GeV/nucleon and CR electrons in the energy range 10 keV- 
10 MeV (hereafter, low-energy CRs) are of particular interest for 
studies of interstellar chemistry, as they provide the bulk of the 
ionization in molecular cloud cores of typical H2 column den- 
sity N ^ 10^^ cm"^. Consequently, low-energy CRs determine 
the degree of coupling of the magnetic field in the gas and con- 
trol the dynamical evolution of molecular cloud cores. Values of 
the ionization fraction in molecular cloud cores determined from 
the column densities of species sensitive to the electron density 
extend over about 1-2 orders of magnitude (Caselli et al. 1998). 
It is unclear whether the inferred CR ionization rates span a 
comparable range or not, due to the sensitivity of the results to 
the adopted chemical model (Williams et al. 1998, Wakelam et 
al. 2006). Since the scatter in the ionization fraction may, in part. 



reflect intrinsic variations of the CR flux from core to core, it is 
important to assess the question of the penetration of CRs in 
molecular cores. 

In a previous study, Padovani et al. (2009, hereafter P09) 
computed the decrease of the CR ionization rate as function of 
column density in a plane-parallel, unmagnetized cloud, consid- 
ering energy losses due to elastic and non-elastic collisions of 
CRs with the ambient gas, bremsstrahlung, and pion production. 
However, interstellar magnetic fields can aflfect the propagation 
of CRs by magnetic focusing, mirroring, and diflTusion (Skilling 
& Strong 1976, Cesarsky & Volk 1978, Chandran 2000, Padoan 
& Scalo 2005). The first two eflTects are due to the non- 
uniformity of the large-scale (mean) component of the field, 
whereas the latter is associated to magnetic field fluctuations on 
the scale of the Larmor radius of CR particles. The relative im- 
portance of these processes in the ISM depends on a number 
of variables not always well determined (geometry and strength 
of the magnetic field, nature and characteristics of turbulence, 
etc.), and is unclear whether they can significantly reduce (or en- 
hance) the CR ionization of molecular cloud cores. In addition, 
the damping of small-scale magnetic fluctuations (e.g. Alfven 
waves) that aflTect the propagation of CRs is strongly depen- 
dent on the ionization fraction of the medium, which, in turn, 
is mostly determined by the CRs themselves. The general prob- 
lem must then be addressed in a self-consistent way. 

In this paper we consider the propagation of CRs along the 
magnetic field lines threading a molecular core assuming that 
only the mean field aflfects the CR properties through a combi- 
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nation of magnetic focusing and mirroring. A similar problem 
was addressed by Desch et al. (2004) to determine the efficiency 
of spallation reactions and trapping of CR nuclei in the proto- 
solar nebula. In their study, Desch et al. (2004) computed an 
energy dependent correction factor for the CR flux, represent- 
ing the fraction of the core's column density a CR of energy E 
interacts with as it passes through the core or it is thermalized. 
Our approach is similar, but we extend and refine their method, 
incorporating the magnetic effects in the propagation model of 
P09, that takes into account energy losses, and adopting a spe- 
cific profile for the density and magnetic field of the core. 

This paper is organized as follows: in Sect. |2]we discuss the 
effects of the magnetic field on the CR propagation; in Sect. [3] 
we describe the adopted cloud core model; in Sect. |4]we com- 
pute the modifications of the CR ionization rate due to the mag- 
netic field; in Sect. [5] we explore the dependence of the results 
on the mass and magnetization of the core; finally, in Sect. [6] we 
summarize our conclusions. 



2. Effects of magnetic fields on CR propagation 

CRs are charged particles that perform helicoidal trajectories 
around the magnetic field lines of the medium where they prop- 
agate. For a particle of charge Ze and mass m in an uniform 
magnetic field, this motion consists of the combination of a cir- 
cular motion around the field line with velocity and a uniform 
motion along the field line with velocity vy. The frequency of the 
circular motion is the cyclotron frequency , Qc = ZeB/ymc. The 
Larmor radius tl is given by 



_pc_ 
ZeB 



sma . 



(1) 



where p = ymv is the particle's momentum and a is the pitch 
angle (i.e. the angle between the particle's velocity and the mag- 
netic field). The helical step S, that is the projection along the 
direction of the magnetic field of the pathlength of the CR dur- 
ing a single rotation around the field lines, is given by 
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For a cloud with a magnetic field B = 10 yuG, the Larmor 
radii of ionizing CRs (CR protons and heavy nuclei with E < 
1 GeV/nucleon and CR electrons with E < 10 MeV, see P09) 
are less than ~ 10"^ pc and ~ 10"^ pc for protons and elec- 
trons, respectively, many orders of magnitude smaller than the 
typical size of Bok globules (~ 0.05 pc), dense cores (~ 1-5 pc), 
and giant molecular clouds (~ 25 pc). In the absence of small- 
scale perturbations in the field, we can therefore assume that 
low-energy CRs propagate closely following the magnetic field 
lines. 

Perturbations in the forms of magnetohydrodynamic (MHD) 
waves with wavelength of the order of the Larmor radius of 
the particle can efficiently scatter CRs. The waves can be part 
of an MHD turbulent cascade, or can be self-generated by the 
CRs themselves (Kulsrud 2005). However, in a mostly neu- 
tral ISM, turbulent MHD cascades are quenched at scales of 
roughly the collision mean free path of ions with neutrals, if 
the ion-neutral collision rate exceeds the energy injection rate. 
Alfven waves with wavelength A < Acr = 7iVA/(yinPi), where 
va is the Alfven speed in the neutrals and 7/„ the collisional 
drag coefficient, are efficiently damped by collisions with neu- 
trals (Kulsrud & Pearce 1969). For typical molecular cloud con- 
ditions (n = 10^ cm"^, B = 10 yuG, and ionization fraction 



~ 10 ^), the critical wavelength is Acr ~ 3 X 10 pc, assuming 
yin ^ 10^"^ cm^ s"^ g"^ (Pinto et al. 2008). As we have seen, the 
Larmor radii of ionizing CR particles are much smaller than the 
cloud's size or the critical wavelength. Thus, for typical values 
of the cloud's parameters, only CR particles with energy larger 
than a few TeV, and Larmor radii tl > Acr find MHD waves to 
resonate with. These particles however do not contribute signif- 
icantly to the ionization of the cloud. 

The origin and maintenance of MHD disturbances at the 
scale of the Larmor radius is unclear: small-scale MHD waves 
can be self-generated by CRs streaming in the intercloud 
medium (ICM) or in the cloud itself. Cesarsky & Volk (1978) 
have shown that the streaming instability operates in the ICM 
only for CRs with energy below ~ 40 MeV/nucleon (for a re- 
gion of size 1 pc, density 10^ cm"^, and magnetic field 10 fiG) 
and at even lower energies in the regions surrounding molecu- 
lar cloud cores. Thus, we expect that CRs in the energy range 
between 100 MeV/nucleon and 1 GeV/nucleon, that provide the 
bulk of the ionization in a cloud core, stream freely through the 
core without self-generating MHD waves. We ignore therefore 
the presence of self-generated waves in cloud cores for the rest 
of this paper, but we believe that this problem deserves further 
scrutiny. 



2.1. Magnetic focusing and mirroring 

Theoretical models predict that collapsing cloud cores must 
overcome the support provided by their magnetic field in order 
to form stars. In the process, the competition between gravity 
pulling inward and magnetic pressure pushing outward is ex- 
pected to produce a warped, hourglass pattern of the magnetic 
field. Recently, this scenario has received support from observa- 
tions. Maps of polarized dust emission have revealed that the 
magnetic field in molecular clouds is rather uniform, except 
near cores where the field becomes strongly pinched and al- 
most radial (see, e.g.. Tang et al. 2009). On the other hand high- 
resolution interferometric observations of submillimeter polar- 
ized emission in the low-mass core NGC 1333 IRAS4A by 
Girart et al. (2006) show a magnetic field geometry consistent 
with the predictions of theoretical models for the formation of 
solar- type stars, in which ordered large-scale magnetic fields 
control the evolution and collapse of molecular cloud cores (see 
a comparison of observations with theoretical collapse models in 
Gonial ves et al. 2008). We therefore adopt the hourglass geom- 
etry as the basis of our analysis of CR penetration into a cloud 
core (see Sect.O. 

The effects of magnetic mirroring and focusing in a hour- 
glass geometry can be simply described following e.g. Desch et 
al. (2004). A charged particle traveling in a magnetized medium 
conserves its kinetic energy ymc^ and its magnetic moment 
= ymv^ sin^ a/2B. It follows that CRs propagating from the 
ICM to the cloud's interior must increase v]_ to conserve jd and 
decrease vy to conserve |vp. Thus, the pitch angle of the particle 
must increase from the value o^icm to a value a as 



sin^ a 



sm a^icM 



B 



5i 



ICM 



(3) 



where ;^ > 1 . Therefore, a CR starting in the ICM with a pitch 
angle sin aicu > 1 /x^^^ cannot penetrate a region with magnetic 
field B > x^iCM, and will be bounced out (magnetic mirror- 
ing). Conversely, the CR flux j(E) in the cloud is increased by 
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Fig. 1. Isodensity contours (thick curves) and magnetic field 
lines (thin curves) in the meridional plane of a model with 
A = 2.66 (Ar = 1.94), from Li & Shu (1996). We consider the 
motion of CRs along a particular field line (red line) as function 
of the polar angle 6. 



the opening out of the field lines by a factor proportional to the 
density of magnetic field lines per unit area (magnetic focusing), 



j(E)=xhcySE). 



(4) 



The eflfects of focusing and mirroring depend only on the mag- 
netic field strength, and are the same for CR protons, electrons, 
and heavy nuclei. 



3. The cloud model 

In order to study the eflfects of the magnetic field on the prop- 
agation of CRs in molecular cores, we adopt the models of 
Li & Shu (1996) (see also Galli et al. 1999), for magneto- 
static, scale-free, self-gravitating clouds supported by axially- 
symmetric hourglass-like magnetic fields (see Figure [T). For 
simplicity, we consider models with an isothermal equation of 
state. These models are characterized by a value of the non- 
dimensional mass-to-flux ratio A defined by 



A = 27iG 



(5) 



where G is the gravitational constant, O the magnetic flux, and 
M(0) the mass contained in the flux tube O. A related quantity 
is the spherical mass-to-flux ratio Ar < A, evaluated from Eq. (O 
but considering only the mass enclosed in a sphere tangent to the 
"waist" of the flux tube. The latter quantity is more appropriate 
for comparisons with observations, because of the limited beam 
of a telescope. 

The magnetic field of the model is axially symmetric and 
purely poloidal. In spherical coordinates (r,6), it is given by 



B 



1 V7 / ^ 

— Vx ( 

In \ r sm ^ 



where 



(D(r, 6) ■■ 



Anc]) 



Cs is the sound speed, and ^(6) is a dimensionless function. 
Similarly, the density is given by 



p(r,^) = 



(8) 



where R(0) is a dimensionless function. From Eq. © and d?]), 
the radial and polar components of the magnetic field are re- 
spectively 



2c 



^0' 



and 



Bb = — 



2c? 



where (p' = AcpjAO. The modulus of the magnetic field vector is 



(9) 



|B| = 



.4j,2 



G3/2m(<D) sin 61 



1 + — + BiCM , 



(10) 



where ^icm is the intercloud magnetic field, assumed uniform 
and equal in strength to the Galactic magnetic field. 

For our reference model we choose A = 2.66 (Ar = 1.94), and 
we assume Cs = 0.2 km s"^ ^icm = 3 fiG. Figure [2] shows the 
profiles of the magnetic field strength and the density as function 
of the polar angle, moving along flux tubes enclosing various 
masses. In the following we will focus on a flux tube enclosing 
a mass M(0) = 1 M©, a typical value for a low-mass core, cor- 
responding to an equatorial radius at the "waist" of the flux tube 
Teq = 0.036 pc (other cases will be considered in Sect.O. 

The upper panel of Figure [3] shows the variation of the pitch 
angle a computed from Eq. (O as function of the polar angle 
6 for various inital values of aicu assuming the magnetic field 
profile of the reference flux tube (see Figure |2). In this case, only 
particles starting with pitch angle q^icm < o^cr = 20.5° are able 
to reach the cloud's midplane (6 = 7r/2, where ;^ = 8.174). CRs 
starting with pitch angles larger than this value will be pushed 
back by magnetic mirroring before reaching the midplane at a 
position ^max(o^iCM) < 7i/2. Inverting this relation, one finds the 
value of the maximum allowed pitch angle CKicM,max(^) for a CR 
to reach a given position 6, shown in the lower panel of Figure [3] 

4. CR ionization rate 

Consider first CRs entering the core from above (in the following 
denoted by a subscript -h). Only those with pitch angle aicu < 

^ICM,max 

(6) can reach a position 6 inside the cloud. Thus, their 
contribution to the CR ionization rate of H2 is 



^^\6) 



27ix(e) dE j[E,N(e,aicM)]^ 
Jo Jo 

X [1 + 0(£)] (t'^'^E) sin aicM daicM , (11) 

where N(6, aicu) is the column density of H2 into the core (with 
A/^ = at ^ = 0), (p(E) is a correction factor accounting for the 
ionization of H2 by secondary electrons, and cr^'^^(E) is the ion- 
ization cross section of H2. We adopt the continuous- slowing- 
down approximation to relate the flux of CR particles of energy 
E 3.13. column density N inside the cloud, j[E, N(0, ^icm)], to the 
flux of CR particles of energy Eiqm in the ICM, j(EicM, 0): 



(6) j[E,N(e,aicM)] = j(EicM,0) 



L(Eicm) 
L(E) ' 



(12) 



where L(E) is the energy loss function (see P09 for details on the 
method). The initial and final energies are such that 



(7) N = n[R(EKM) - R(E)], 



(13) 
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Fig. 2. Modulus of the magnetic field (upper panel) and density 
(lower panel) as function of the polar angle 6 for flux tubes en- 
closing different masses (label values in M©). 



where n is the local density of H2 and R(E) is the particle's 
range. We consider CRs composed of electrons, protons, and 
heavy nuclei. As stressed by Webber (1998) and P09, CR elec- 
trons contribute significantly to the total ionization rate, and can- 
not be neglected. We adopt four combinations of proton/electron 
spectra j(E), labeled as in P09: Webber (1998, W98) and 
Moskalenko et al. (2002, M02) for protons and heavy nuclei; 
Strong et al. (2000, COO and EOO) for electrons. The spectra 
M02 and EOO are characterized by an increase at low energies, 
whereas the spectra W98 and COO are constant or slightly de- 
creasing. 

The column density passed through by a CR propagating 
along a magnetic field line is 



N(e, aicu) 



cos a 



(14) 



where [i = 2.8 is the molecular weight, d£ = [dr^ -h (rd6)^Y^^ 
is the element of magnetic field line, and the factor 1 / cos a ac- 
counts for the increase of the actual path length of a CR perform- 
ing a helicoidal trajectory around a magnetic field with respect to 
the displacement along the field line. Substituting in Eq. ([T4l) the 
expressions for p and d£, of our model, we obtain the actual col- 
umn density passed through by CRs propagating from the ICM 
to the mirror point ^max(a^iCM), 



A^(^, aicu) = 



7i/dmnG^M((^) 



re, 
Jo 



mm 

cos a(6) 




0.4 0.6 
e [n/2 rad] 

Fig. 3. Upper panel: variation of the CR pitch angle a as function 
of the polar angle 6 for values of o^icm between 0° and 90° in 
steps of 5° (from bottom to top); the dashed line represents the 
critical pitch angle a^r = 20.5°. Lower panel: relation between 
the maximum initial pitch angle aicM,max and the polar angle that 
can be reached during the propagation inside the core. 



where cos a(6) is given in terms of OficM and;^ by Eq. (O. 

The results are shown in Figure |4] for the reference flux tube 
enclosing 1 M© and for different values of aicu- It is evident 
from the Figure that as a first approximation one can assume 
that CRs coming from the ICM and traveling toward the cloud's 
midplane experience a similar increase in column density, inde- 
pendently on the initial pitch angle aicu, the latter mainly de- 
termining the value of the column density at which the CRs are 
pushed out by magnetic mirroring. We will then assume for all 
CRs the column density profile corresponding to o^icm = 0, but 
truncated at increasingly larger values depending on the initial 
pitch angle OfjcM, 



A^(^,^^icm)^A^(^,0) if 0<^<^,,ax(^:^icM). 
With this approximation Eq. (fTTI) simplifies to 



(16) 



f 

Jo 



ric 



.x(^) 



sin aicM dcKicM x 



j[E,No(e)] [l+0(£)]cr^^^(£) dE, 



(17) 



+ P- d^,(15) 



where, for simplicity, we have defined A^o(^) = N(6, 0) to indi- 
cate the column density measured along a field line. The second 
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Fig. 4. Column density passed through before reaching the mir- 
ror point by CRs propagating along a field line of the reference 
flux tube enclosing 1 M©, as function of the polar angle 6 for dif- 
ferent values of the initial pitch angle OficM (labels in degrees). 



integral of Eq. (fTTl) is the CR ionization rate computed assuming 
one-dimensional propagation, calculated in P09: 



Jo 



A^o(^)] [1 + (pm (T^'^'XE) dE ^ —^^^[Nom . 



471^ 



Thus, Eq. (fTTl) can be rewritten as 
where 



(f+(0) = - COSaiCM,n 



(18) 



(19) 



(20) 



is the factor accounting for magnetic eflTects. This completes the 
calculation of the contribution to the ionization of CRs coming 
from the upper side of the cloud. 

A CR propagating along a field line and reaching a polar an- 
gle 6 coming from the lower side of the cloud has passed through 
a column density 



No(7t/2) + [No(7t/2) - Nom = 2No(7t/2) - No(6) , 



(21) 



where No(7:/2) is the cloud's column density at the midplane. 
Thus, the contribution of CRs coming from the lower side of the 
cloud (denoted by a subscript -) to the ionization rate is 



where ip-{6) is 

^_(0=^(1- cosher). 

The total CR ionization rate is then given by 



(22) 



(23) 



(24) 



Figure [3 shows ^^^{6) including the eflTects of magnetic focus- 
ing and mirroring, compared to the CR ionization rate evaluated 
without considering the magnetic field. The curves refer to the 
four combinations of the two protons and electron spectra as- 
sumed in this work (see the paragraph following Eq. [T3]l. It is 



evident that the two processes have opposite eflTects of compa- 
rable magnitude on the ionization rate, both becoming more and 
more important approaching the core's midplane where the field 
is stronger. In Figure |6] we compare the values of computed 
with and without magnetic eflTects as function of the column den- 
sity (measured along field lines) in the core. The quantity % de- 
fined as the ratio of the rates obtained in the two cases, is shown 
in the lower panel of the Figure. Clearly, magnetic mirroring al- 
ways reduces the CR ionization rate more than magnetic focus- 
ing can increase it, the total eflfect being a net reduction of f 
by a factor between 2 and 3. Notice that the maximum eflTect 
of the magnetic on ^gt obtained on the core's midplane, 

where the field is stronger, but at an intermediate position cor- 
responding to column densities 10^^-10^^ cm"^. In fact, in the 
core's midplane ^+(7r/2) = ^_(7r/2) = ;^(7r/2)(l - cosacr)/2, 
and, using Eq. (O, one has 



n{7ll2) =xi7^m{^ - cosher) = 1^) 



1- Jl- 



1 



Xi7il2) 



(25) 



independent on the assumed CR spectrum, a result also obtained 
by Desch et al. (2004). For the reference flux tube enclosing 
1 Mo, ^(7r/2) = 0.516. 

5. Dependence of on the mass and 
magnetization of tlie core 

5.1. Varying the enclosed mass M(0) 

In the previous section we have considered CR propagating 
along the field lines of a flux tube enclosing 1 M©. In Table [T] 
we show the variation of the relevant quantities of the prob- 
lem as function of the enclosed mass for our reference core with 
A = 2.66. Flux tubes enclosing smaller masses intersect the mid- 
plane at smaller req and are characterized by larger values of 
the magnetic field and density (see Figure In particular, in 
our model the field strength and the density increase inward as 
|B| oc r"^ oc M(0)"^ and p oc oc M(0)"^, respectively (see 
Eq. [TOland Eq. (8]). As a consequence, focusing becomes more 
important, since x increases, but also mirroring becomes more 
severe, since o^cr decreases. The net effect is a stronger reduc- 
tion of the innermost regions of the core as compared to 
the envelope (see Table [TJ. As the field strength increases ap- 
proaching the central singularity of the model, x ^ ^ and the 
reduction factor of the CR ionization rate in the core's midplane 
approaches the asymptotic value 1i{7il2) =1/2. 

Conversely, for flux tubes enclosing larger masses, the field 
strength approaches the ICM value, and the density decreases 
to zero. Therefore both focusing and mirroring become weaker, 
and o^cr approaches 90°, as shown by Table [T] As expected, for 
increasing values of M(0), approches the value of the non- 
magnetic case, because the magnetic field strength decreases 
away from the center of the core, approaching the ICM value. 

5.2. Varying the mass-to-flux ratio A 

In this section we explore the effects of the variation of the 
core's mass-to-flux ratio A on for a flux tube containing 
a fixed mass of 1 M©. For cores with strong magnetic support 
(lower values of A), the equatorial squeezing of the field lines is 
stronger, and the lines reach more internal regions of the core 
where the density is higher. As a consequence, acr decreases, the 
mirroring effect becomes stronger, and a smaller fraction of CRs 
can penetrate the cloud. In Figure |7] we show the ratio as a 
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0.4 0.6 
Q [n/2 rad] 

Fig. 5. Effects of magnetic mirroring (upper panel, solid lines) 
and focusing (lower panel, solid lines) on the CR ionization 
rate as function of the polar angle. Dotted lines represent the 
case without considering the magnetic field influence, from P09. 
Labels refer to the CR proton and electron spectra used in P09 
in the non-magnetic field case (proton spectra: Webber 1998, 
W98; Moskalenko et al. 2002, M02; electron spectra: Strong et 
al. 2000, COO and EOO). 



Table 1. Values of the parameters described in the text as func- 
tion of the mass, M(0), contained within flux tube O. 



M(0) 




No(n/2) 




X(7r/2) 


m/2) 


(Mo) 


(pc) 


(10^3 cm-2) 








0.5 


0.018 


4.54 


14.8° 


15.347 


0.508 


1 


0.036 


2.27 


20.5° 


8.174 


0.516 


5 


0.180 


0.45 


39.9° 


2.435 


0.566 


10 


0.360 


0.23 


49.7° 


1.717 


0.608 


50 


1.802 


0.05 


69.3° 


1.143 


0.739 


100 


3.604 


0.02 


75.0° 


1.072 


0.794 



function of column density for decreasing values of A. Eor sim- 
plicity, we have considered only the combination of the proton 
spectrum M02 and the electron spectrum COO (the results ob- 
tained with the other spectra are similar). As the Eigure shows, 
the reduction of larger in cores with larger magnetic sup- 
port, due to the increase in the field strength and concentration 
of field lines. The reduction is a factor ~ 4 for the outer regions 
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Fig. 6. Comparison between the CR ionization rate with and 
without the effects of magnetic field (solid and dotted lines, re- 
spectively) as function of the column density (upper panel)', ratio 
between the ionization rates in the magnetic and non-magnetic 
case (lower panel). The curves are labeled as in Eigure [5] 



of cores with A = 1.63, the lowest value of mass-to-flux ratio 
considered in our models. 

Eor A ^ I, the density distribution becomes more and more 
flattened, the core assumes the shape of a thin disk, and the 
column density from the ICM to the core's midplane becomes 
larger. In this limiting case, the magnetic field strength and the 
column density increase as (A^ - 1)"^. Eor these magnetically 
dominated disk- like configurations, the reduction of the CR ion- 
ization rate approaches the asymptotic value 1/2. 

6. Conclusions 

We have determined the CR ionization rate f a magneti- 
cally supported molecular core, and compared the results with 
the case of non-magnetized density distribution with the same 
column density. We have considered only the large-scale mag- 
netic field threading the core because, according to previous in- 
vestigations, CR scattering by magnetic fluctuation on the scale 
of the Larmor radius of the particle should not be important for 
CRs in the energy range 100 MeV-1 GeV, responsible for the 
bulk of the ionization and heating in molecular clouds. 

Assuming the hourglass magnetic field profile and the den- 
sity distribution of a model of magnetically supported cloud 
core, we computed function of position in the core, for dif- 
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Fig. 7. Ratio between the CR ionization rates in the magnetic 
and non-magnetic case for the case of M02-hC00 spectrum (see 
Figure [3 for the labels). The curves are computed for a flux tube 
containing 1 M© and for difl'erent values of A = 8.38, 2.66, 
2.07, 1.79, 1.63. The upper scale shows the extinction through 
the cloud obtained from Ay /mag = A^/10^^ cm"^. 



ferent values of the intercloud CR flux and composition, taking 
into account the energy loss processes considered in the previous 
work of P09 and following the CR propagation along magnetic 
field lines. The relatively weak dependence of the reduction of 
the CR ionization rate on the mass enclosed in a flux tube (basi- 
cally a logarithmic dependence, see Tab.[T]), and on the mass-to- 
flux ratio of the cloud core (A ^ 2on average for most cores, see 
Troland & Crutcher 2008), lends some confidence on the general 
validity of our results beyond the specific cloud model adopted 
here. The method described in this paper can be applied however 
to any axisymmetric cloud model, provided the input quantities 
B and p along a given flux tube are specified as function of the 
polar angle or vertical cylindrical coordinate, as in Fig. |2] 

Our results show that f a magnetized core is always re- 
duced with respect to its non-magnetic value, by a factor de- 
pending on the core's mass-to-flux ratio (determining the field 
strength) and the amount of mass contained in the flux tube con- 
sidered. In general, for a flux tube enclosing about 1 M©, and for 
a core with mass-to-flux ratio A = 2.66, the CR ionization rate is 
reduced by a factor of ~ 3 over most of the core, and by a factor 
of ~ 2 near the core's midplane. The reduction is less severe for 
flux tubes enclosing larger masses and for larger values of the 
core's mass-to-flux ratio, because of the near uniformity of the 
magnetic field in these cases. For cores with mass-to-flux ratio 
A ^ 1.6 or lower, the magnetic field can decrease the CR ioniza- 
tion rate by a factor of ~ 4 over most of the core's envelope. 

Thus, the values of derived by Caselli et al. (1998), 
WilHams et al. (1998), Maret & Bergin (2007) for dense cores 
and globules, probably underestimate the "external" (i.e. inter- 
cloud) CR ionization rate by a factor of ~ 3 - 4, thus alleviating 
the discrepancy with measurements of diflfuse clouds. In 

this respect, it is interesting to notice that the best overall fit with 
the chemical abundances of the low-mass core TMC-1 is ob- 
tained when the CR ionization rate is reduced by a factor ~ 5 
with respect to the standard value ~ 10"^^ s"\ together 
with the rate of removal of atomic hydrogen from the gas phase 
(Flower et al. 2007). However, the observed large scatter in the 



values of of dense clouds found by Caselli et al. (1998), if 
real, cannot be attributed to eflfects of the mean field or to varia- 
tions in the mass-to-flux ratio of the cores. 
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